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Summary  of  Progress  and  Results 

1.  Photodissociation  dynamics  of  CIN3 

Photolysis  of  C1N3  is  a  useful  source  of  NCl(a)  for  kinetic  measurements.  For  experiments 
conducted  under  pseudo  fist  order  conditions  it  is  sufficient  to  know  the  relative  concentration  of 
NCl(a),  but  the  absolute  value  is  not  needed.  However,  to  determine  rate  constants  under  second 
order  conditions  knowledge  of  the  absolute  concentration  of  NCl(a)  is  needed.  One  of  the  goals  of 
this  research  has  been  measurement  of  the  second  order  rate  constant  for  self-annihilation  of  NCl(a). 
The  technique  we  have  applied  to  generate  known  quantities  of  NCl(a)  involves  complete 
photodissociation  of  a  known  concentration  of  C'1N3.  This  can  yield  well-defined  quantities  of  NCl(a) 
provided  that  the  branching  fraction  for  this  product  is  known.  It  is  established  that  NCl(o)  is  the 
dominant  product,  but  NCl(fe)  and  NC1(X)  are  also  produced.  In  addition,  the  channel 
CINj+hv— >C1+N3  is  open.  Previously  we  had  used  dispersed  fluorescence  to  observe  singlet  NCI  and 
absorption  measurements  to  monitor  NC1(X).  The  limitations  of  this  approach  were  that  we  could  not 
verify  total  destruction  of  C1N3  by  direct  observation,  and  that  non-fluorescing  products  (other  than 
NC1(X))  could  not  be  monitored.  To  improve  our  knowledge  of  C1N3  photodissociation  dynamics  we 
examined  the  potential  for  using  photo-ionization  mass  spectrometry  as  a  more  versatile  detection 
method.  We  began  these  studies  by  using  pulsed  118  nm  light  to  ionize  C1N3.  To  generate  VUV 
pulses  we  added  a  frequency-tripling  cell  to  our  TOF  mass  spectrometer,  which  was  used  to  up 
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convert  the  third  harmonic  from  an  Nd/YAG  laser.  Fig.  1  shows  the  mass  spectrum  obtained  using 
118  nm  pulses  to  excite  C1N3.  Surprisingly,  the  azide  did  not  survive  this  rather  gentle  method  of 
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Figure  1.  Mass  spectrum  obtained  by  1 18  nm  excitation  of  C1N3 


ionization.  The  position  where  the  C1N3  should  appear  is  indicated  in  Fig.  1,  but  there  was  no 
discemable  signal  at  this  mass.  It  appears  that  118  nm  light  causes  a  dissociative  ionization  of  C1N3 
of  the  type 

CINj  +  hv(l  18  nm)  NCT  +  N2. 

Hence  this  approach  is  not  suitable  for  monitoring  C1N3  in  experiments  where  alternative  pathways  for 
generating  NCI  exist.  We  also  examined  the  possibility  of  using  excimer  laser  pulses  (193  or  248  nm) 
to  effect  multi-photon  ionization  of  C1N3  The  results  of  these  experiments  are  shown  in  Figs  2  and  3. 
Note  that  the  excimer  laser  pulses  fragmented  and  ionized  the  background  pump  oil  in  the  mass 
spectrometer  chamber  (evacuated  to  2xl0'7  Torr).  The  downward-going  traces  of  Figs  2  and  3  show 
the  background  spectra,  which  consisted  mostly  of  hydrocarbon  fragments.  The  upward-going  traces 
were  taken  with  C1N3  present.  Neither  of  the  photoexcitation  wavelengths  generated  detectable 
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Figure  2.  Mass  spectrum  obtained  by  multiphoton  excitation  of  C1N3  at  193  nm 


Figure  3.  Mass  spectrum  obtained  by  multiphoton  excitation  of  C1N3  at  248  nm 


quantities  of  C1N3\  Again,  fragmentation  to  produce  NC1+  dominated  the  kinetics.  Note  that 
excitation  at  193  nm  produced  N3+,  which  suggests  the  possible  production  sequence 

C1N3  +  hv  — )  Cl  +  N3 
N3  +  2hv  N3+ 

Evidence  that  193  nm  light  can  cause  the  first  step  was  obtained  in  a  study  of  the  IR  emission  (see 
below).  The  arrow  in  Fig.  3  marks  the  position  where  N3  should  be  observed  in  the  mass  spectrum 
recorded  using  248  nm  light.  The  absence  of  N3  in  this  spectrum  may  indicate  that  the  path  to  Cl+N3 
represents  a  very  minor  channel  for  248  nm  photolysis. 

In  previous  studies  we  examined  the  kinetics  of  NC1(X)  production  when  ClN3/02  mixtures 
were  photolyzed.  The  purpose  of  these  measurements  was  to  determine  the  branching  fraction  for 
fonnation  of  NCl(a).  As  the  experiments  involved  monitoring  a  single  ro-vibrational  level  of  NC1(X) 
v=0,  the  role  of  ro-vibrational  relaxation  in  the  observed  kinetics  was  an  open  question.  Initial 
kinetic  models  required  the  assumption  that  UV  photolysis  of  C1N3  yielded  NC1(X)  v>l  with  a 
branching  fraction  of  about  20%.  However,  attempts  to  detect  population  in  the  v=l  and  2  states  by 
transient  absorption  were  unsuccessful.  As  an  alternative  means  for  characterizing  these  products  we 
attempted  to  detect  IR  fluorescence  from  NCI  following  photolysis  of  C1N3.  These  experiments  were 
conducted  in  collaboration  with  Prof.  Hai-Lung  Dai  at  the  University  of  Pennsylvania.  The  apparatus 
used  is  shown  in  Fig.  4.  C1N3  was  photolyzed  by  193  nm  pulses  from  an  excimer  laser.  The 
frequency  and  temporal  characteristics  of  the  resulting  IR  fluorescence  were  recorded  using  a  time- 
resolved  FTIR  spectrometer.  Fig.  5  shows  a  typical  emission  spectrum.  The  only  significant  feature 
in  this  trace  is  the  emission  from  N3  (asymmetric  stretch  fundamental)  at  1644  cm'1.  There  was  no 
sign  of  IR  emission  from  NC1(X).  Unfortunately,  the  absence  of  NC1(X)  bands  is  not  conclusive 
evidence  that  this  product  is  vibrationally  cold.  The  transition  dipole  for  NC1(X)  IR  emission  is  rather 
small,  and  occurs  near  the  long  wavelength  cut-off  of  the  HgCdTe  detector.  Conversely,  the  transition 
dipole  for  the  asymmetric  stretch  of  N3  is  relatively  large,  and  the  emission  occurs  in  a  spectral  range 
that  is  better  suited  for  the  detector.  Hence,  although  we  could  observe  N3this  does  not  necessarily 
indicate  that  dissociation  to  give  CI+N3  is  a  dominant  channel.  By  comparison  with  isovalent  HN3  we 
estimate  that  about  15%  of  the  193  nm  photolysis  events  give  Cl+N3. 
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Figure  4.  Apparatus  used  for  time-resolved  FTIR  study  of  C1N3  photolysis 
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Figure  5.  IR  emission  from  photolysis  of  C1N3  at  193  nm. 
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Different  electronic  transitions  of  C1N3  are  accessed  by  193  and  248  nm  excitation.  A  similar 
situation  exists  for  HN3,  where  the  maxima  of  the  two  systems  are  shifted  to  longer  wavelengths. 
Photolysis  of  HN3  via  the  lower  energy  state  (e.g.,  266  nm  excitation)  gives  a  small  yield  of  H+N3 
(4%),  whereas  excitation  to  the  higher  energy  state  (248  nm  excitation)  yields  20%  H+N3.  The  data 
from  our  mass  spectrometry  experiments  suggest  that  C1N3  behaves  similarly,  and  that  the  yield  of 
Cl+N3  from  248  nm  photolysis  (excitation  to  the  lower  energy  state)  will  be  small. 

In  future  investigations  of  the  IR  emissions  we  plan  to  examine  248  nm  photolysis  and  use 
detectors  that  are  better  suited  for  800-1000  cm4  radiation. 

Combining  a  conservative  estimate  of  15%  for  the  248  nm  yield  of  Cl+N3  with  the  NCl(a) 
branching  fraction  from  the  NC1(X)  absorption  measurements  ([NCl(a)]o/([NCI(a)]0+[NCI(X)]0)>0.9) 
we  obtain  a  quantum  yield  for  formation  of  NCl(a)  of  76%.  This  value,  applied  to  the  data  from  our 
previous  NCl(a)  self-removal  study,  yields  a  second  order  rate  constant  of  6.2xl043  cm3  s4.  Henshaw 
et  al.  reported  a  self-annihilation  rate  constant  that  was  more  than  an  order  of  magnitude  greater  than 
this  value  ((7.2±G.9)xl042  cm3  s4).  For  our  fluorescence  decay  data  to  be  consistent  with  this  value 
the  quantum  yield  for  NCl(a)  formation  would  have  to  be  just  7%.  This  is  unreasonably  low.  All 
other  estimates  for  the  quantum  yield  exceed  50%.  The  demonstration  of  I*  lasing  via  photolysis  of 
C1N3/CH2I2  mixtures  further  attests  to  the  high  quantum  yield. 

2,  Construction  of  a  heated  photolysis  cell 

In  the  original  proposal  we  had  planned  to  study  the  temperature  dependencies  of  the  NCl(a) 
quenching  rate  constants  and  the  self-annihilation  reaction  using  a  heated  flow  tube  with  mass 
spectrometric  detection.  We  have  since  discovered  complications  with  this  approach  that  are 
associated  with  the  dissociative  photoionization  of  C1N3,  and  developed  a  technique  for  characterizing 
die  self-annihilation  process  using  fluorescence  detection.  Consequently,  we  are  now  constructing  a 
heated  photolysis  cell  that  is  configured  for  fluorescence  detection.  This  apparatus  will  be  used  to 
explore  the  temperature  dependence  of  NCI(a)+M,  where  M=NCl(a),  HF,  HC1,  HI,  Cl2,  and  C1N3. 
Recent  model  calculations  (P.  Crowel,  Logicon  RDA)  indicate  that  the  Cl  atoms  released  by  the 
reaction 

NC1(X)  +  NC1(X)  N2  +  2C1  (1) 
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may  pose  a  serious  problem  in  the  NCI/I  laser  as  Cl  is  an  efficient  quencher  of  I*.  At  present  the 
temperature  dependence  of  the  I*+C1  quenching  rate  constant  is  not  known.  The  heated  photolysis 
cell  will  be  used  to  obtain  this  information. 

3.  Theoretical  studies  of  NCI  +  NCI  interactions. 

To  gain  deeper  insights  concerning  the  self-annihilation  of  NCl(a)  and  the  range  of  possible 
products  from  NC1(X)  +  NC1(X)  we  have  examined  NCI  -  NCI  interactions  using  high-level 
theoretical  methods.  This  work  is  being  carried  out  in  collaboration  with  Prof.  Keiji  Morokuma  and 
Dr.  Gregory  S.  Tschumper  at  Emory. 

One  of  the  issues  explored  in  this  study  was  the  stability  of  ground  state  N2C12.  Recent  BP86 
density  functional  calculations  had  indicated  that  both  the  cis-  and  trans-  isomers  were  stable,  with  the 
former  being  a  little  more  deeply  bound.  We  applied  a  number  of  multi-reference  techniques  in  our 
calculations.  The  active  spaces  used  for  these  calculations  can  be  defined  in  terms  of  the  natural 
orbitals  depicted  in  Fig.  6  (for  the  trans-  isomer).  The  6  electron  6  orbital  (6x6),  and  8  electron,  8 
orbital  (8x8)  spaces  are  indicated.  A  few  full  valence  calculations  (24x16)  were  also  carried  out.  The 
6-311+G(d)  basis  set  was  used  with  the  CASSCF  and  CAPT2  methods.  In  addition,  coupled  cluster 
calculations  (CCSD(T))  were  performed  with  the  TZ2P  and  TZ2P (f)  basis  sets.  Results  from  these 
calculations  for  the  cis-  and  trans- isomers  are  presented  in  Fig.  7  and  Table  1.  Here  we  use  the 
notation  N2Cl2f  to  designate  the  saddle  point  that  separates  N2C12  from  N2+2C1.  Regardless  of  the  size 
of  the  active  space,  CASSCF  predicts  that  N2C12  is  unstable  with  respect  to  dissociation  to  N2+2C1.  It 
appears  that  the  previous  predictions  of  stable  structures  were  an  artifact  resulting  from  the  use  of  a 
single  reference  approximation.  CASPT2  and  CCSD(T)  calculations  predict  very  shallow  minima  for 
the  cis-  and  trans-  structures.  The  exact  locations  of  the  CASPT2  saddle  points  were  not  determined. 
The  CCSD(T)/TZ2P(/)  saddle  points  yield  barriers  to  dissociation  of  8.0  kcal/mol  (trans)  and  7.3  kcal 
/mol  (cis).  These  barriers  are  not  high  enough  to  make  N2C12  a  stable  species.  This  is  consistent  with 
the  fact  that  there  have  been  no  experimental  observations  of  this  molecule. 
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Figure  6,  Natural  orbitals  for  trans-  C1N-NC1 
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Energies  (in  kcal/mol)  of  key  stationary  points  relative  to  the  NCI  (.X^E  )  +  NCI  (X^E  )  assymptote. 


Method 

Basis 

imm- N2CI2 

CM-N2G2  IranJi-NaCV 

cifl-N2CV 

n2+  2  a 

6x6  CASSCF 

6-311 +G(d) 

minima  and  transition  states  do  not  exist 

-116.2 

8x8  CASSCF 

6-311+G(d) 

minima 

and  transition  state?  do  not  exist 

*  *  * 

24  x  16  CASSCF 

6411+G(d) 

minima  and  transition  states  do  not  exist 

-137.9 

6x6  CASPT2 

6-311+G(d) 

-82.4 

-84.4 

*  *  * 

♦  ft 

8x8  CASPT2 

6-311 +G(d) 

-83.2 

-86.5 

*  *  * 

+  *  ♦ 

CCSD(T) 

TZ2P 

-86.5 

-92.9  -82.0 

-87.0 

-106.1 

CCSD(T) 

TZ2P  (/) 

-88.6 

-92.6  -80.6 

-85.3 

-102.3 

/ 

Figure  7 


Considering  only  the  first  three  bound  electronic  states  of  NCI  (X3Z\  a1  A  and  fe*Z+),  the 
interaction  of  two  NCI  molecules  yields  18  electronic  states,  A  one-dimensional  cut  through  the  C1N- 
NC1  potential  energy  surfaces  of  the  lower  energy  states  is  shown  in  Fig,  8,  This  plot  corresponds  to 
a  rigid  C2h  geometry,  with  variation  of  just  the  N-N  distance.  An  interesting  feature  of  Fig,  8  is  the 
crossings  between  the  lowest  states  of  the  a*A+a* A  with  states  from  X3E+blE  that  occur  near  2,4  A, 
These  crossings  persist  for  many  other  geometries,  and  they  occur  at  energies  (relative  to  the 
alA+alA  asymptote)  that  are  accessible  to  room  temperature  collisions.  Table  2  lists  the  energies  and 
geometries  of  the  curve  crossing  for  the  cis-  and  trans-isomers. 
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Rigid  C2h  Scan  oftheCIN...NCI  Surfaces 
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Table  1 


Optimized  geometrical  parameters  (in  A  and  degrees)  of  key 
stationary  points  on  the  NO  (X3S“)  4  NCl  (X3E-)  surface. 


Method _ Basis  R(NN)  R(NQ)  6(NNQ) 

Isolated  Ns  and  NCl 


6  x  6  CASFT2 

6-31l40(d) 

f  r  t 

1.626 

8x8  CASPT2 

6-31140(d) 

1.639 

#  9  W 

E/ROOOSD(T) 

TZ2P 

1.103 

1.665 

*  *  * 

e/rooosdct) 

TZ2F  (f) 

1.104 

1.640 

ci?NaOa 

t  r  r 

6  X0OASPT2 

6-31 140(d)  ' 

1.206 

1.800 

122.6 

8x8  CASPT2 

6-31 140(d) 

1.203 

1.803 

122.6 

EOOSDCI) 

TZ2P 

1.201 

1.823 

121.9 

ECOSDCI) 

TZ2F(f) 

1.200 

1.793 

transNsCls 

121.4 

6x6  GASPT2 

6-3H4G(d) 

1.223 

1.802 

108.4 

8  x  8  CASPT2 

6-31 140  (d) 

1.224 

1.807 

103.2 

EOOSDCI) 

TZ2F 

1.221 

1.819 

107.4 

EOCSD{T) 

TZ2P  (f) 

1.227 

1.786 

107.8 

6  X6CASPT2 

6-31140(d) 

r  f  f 

f  f  t 

w  p  w 

8x8  GASPT2 

6-31 140(d) 

9  W  W 

f  f  f 

WWW 

ECOSD(T} 

TZ2P 

1.142 

2.082 

124.3 

EOOSD{T) 

TZ2F(f) 

1.137 

2.115 

frun1,rN2CM 

125.1 

6x6  0ASFT2 

6-31 140(d) 

... 

... 

... 

8  x8  CASPT2 

6-31 140(d) 

... 

... 

... 

EOOSD(I) 

TZ2P 

1.157 

2.086 

105.2 

EOOSD(T) 

TZ2P(f) 

1.150 

2.115 

105.8 
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Table  2 


Optimized  geometrical  parameters  (in  A  and  degrees)  and  rela¬ 
tive  energies  (in  kcal  mol-1 /mol)  of  Rome  long  range  ZUiSXs  on 
tiie  o  +  o  surface.  The  energies  of  the  MSXh  are  relative  to 
the  NO  (a1  A)  +  NO  (oxA)  assyraptotes.  All  value:  rare  ob¬ 
tained  using  the  16  X  12  state -averaged  conjunctions  with  the 
6-311+Q(d)  basis  set. 


X  +  6 

a  +  o 

R(NN) 

R(NO]) 

0(NNC1) 

A-Eat5X 

d,?-Oonjbimstion 

T7 

SI  (M3) 

2.373 

1.634 

114.3 

3.3 

T7  0BX) 

S2  (%) 

2.351 

1,632 

114,4 

4.3 

h^ri^Oonfcormation 

T 7<%) 

S1W 

2.345 

1.643 

10942 

3.2 

W(M4 

saw 

2,375 

1.645 

109.7 

2.2 

It  is  probable  that  these  crossings  are  responsible  for  the  energy  pooling  process 
NCl(a)+NCl(a)  ->  NCl(b)+NCl(X) 

Based  on  our  NCl(a)  self-annihilation  rate  constant  and  the  rate  constant  for  reaction  2  reported  by 
Manke  and  Setser,  the  pooling  reaction  accounts  for  approximately  25%  of  the  self  removal  process. 
The  calculations  performed  so  far  support  the  notion  that  reaction  1  rapidly  removes  NC1(X) 
and  they  indicate  that  electronic  energy  transfer  via  curve  crossings  contributes  significantly  to  the 
self-removal  of  NCl(a),  However,  further  test  are  needed  to  ensure  that  we  have  adequately 
converged  results.  Selected  calculations  that  use  full-valence,  state  averaged  MCSCF  will  be  used  to 
test  the  reliability  of  the  lower  level  methods. 
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5.  Theoretical  studies  of  the  Cl+N3  —>  NCI  +  N2  reaction. 

The  reaction  Cl  +  N3  -¥  NCI  +  N2  is  known  to  give  substantial  yields  of  NCl(a).  The  potential 
use  of  this  reaction  in  a  chemically  driven  iodine  laser  system  is  currently  being  evaluated.  The 
branching  ratio  for  Cl  +  N3  — »  N2  +  NCl(a,  b)  (singlet  product)  vs.  Cl  +  N3  — »  N2  +  NC1(X)  (triplet 
product)  is  a  critical  issue  for  laser  pumping  applications.  In  recent  kinetic  studies  it  has  been 
determined  that  NCl(a)  is  formed  with  a  yield  in  excess  of  60%.  However,  it  is  a  difficult  matter  to 
make  a  direct  experimental  determination  of  the  branching  ratio.  To  address  this  issue  we  are 
performing  ab  initio  calculations  to  examine  the  singlet/triplet  branching  ratio  from  a  theoretical 
perspective.  This  work  is  being  done  in  collaboration  with  Prof.  Keiji  Morokuma  and  Mr.  Graham 
Phillips. 

Previously  we  examined  qualitative  features  of  the  Cl-N3  potential  energy  surfaces  using  the 
B3LYP/6-3 1+G(<2)  level  of  theory.  Improved  energies  for  the  stationary  points  were  then  computed 
using  CASSCF  (12  electrons,  10  orbitals)  followed  by  CASPT2  (8  electrons,  6  orbitals).  The  basis  set 
used  was  of  valence  triple  zeta  quality  (spdf,  sp+).  A  puzzling  feature  of  the  results  from  these 
calculations  was  the  low  barrier  (below  the  asymptotic  reactant  energy)  for  the  reaction  on  the  triplet 
surface  (Cl+N3  — >  NC1(X)  +  N2).  With  such  a  low  barrier  it  was  difficult  to  account  for  the  observed 
preference  for  the  formation  of  singlet  products.  We  have  now  investigated  this  system  using  better 
quality  basis  functions.  Our  ability  to  perform  such  large  calculations  was  greatly  improved  by 
upgrades  in  the  Emerson  Center  computing  facilities  that  were  implemented  at  the  beginning  of  this 
year.  New  maps  of  potential  surfaces  were  constructed  using  B3LYP  with  the  D95  +d  basis. 
Stationary  points  and  minimum  energy  crossing  points  were  optimized  using  CASSCF  (12e/10o)  with 
the  D95  +d  basis.  Results  from  the  CASSCF  calculations  are  summarized  in  Fig.  9.  They  are 
dramatically  different  from  the  results  obtained  with  the  6-31+G(d)  basis.  We  now  find  that  the 
barrier  to  reaction  on  the  triplet  surface  lies  appreciably  above  the  entrance  channel,  while  the  barrier 
to  reaction  on  the  singlet  surface  lies  below  the  energy  of  the  reactants.  This  readily  explains  the 
observed  preference  for  the  singlet  products.  Note  that  a  seam  of  intersection  exists  between  the 
singlet  and  triplet  surfaces,  and  the  minimum  of  this  seam  is  indicated  by  the  point  S0T1TX  in  Fig.  9. 
This  intersection  occurs  before  the  barrier  on  the  singlet  surface,  and  may  well  influence  the 
singlet/triplet  product  branching  ratios. 
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Next  we  will  carry  out  dynamical  studies  using  the  potential  energy  surfaces  to  determine  the 
product  branching  ratios.  Calculations  are  underway  to  compare  the  reaction  probabilities  on  the 
singlet  and  triplet  surfaces.  These  are  classical  calculations  that  do  not  include  the  influence  of  the 
intersection  seam.  In  subsequent  calculations  the  coupling  between  the  triplet  and  singlet  states  will 
be  evaluated,  and  the  dynamics  investigated  using  surface  hopping  trajectory  methods. 


Energy, 

kcal/moi  Figure  9:  CASSCF(12e/10o)/D95+d  Potential  Energy  Surface 
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6.  Measurement  of  the  rate  constants  for  H+F2  and  H+Cl2 

Investigation  of  the  rate  constants  for  H+F2  and  H+C12  was  initiated  by  Dr.  G.  Manke  II. 
(AFRL/DELC).  After  performing  an  extensive  review  of  the  kinetic  data  for  the  HF  chemical  laser, 
Dr.  Manke  found  that  the  uncertainty  associated  with  the  rate  constant  for  the  reaction  H+F2— >HF+F 
was  unacceptably  large  (models  of  the  laser  system  are  strongly  dependent  on  this  parameter).  As  the 
facilities  at  Emory  are  better  suited  than  those  at  AFRL/DELC  for  measurements  of  this  kind,  we  have 
collaborated  on  new  studies  of  H  atom  kinetics.  The  H+C12  reaction  was  investigated  as  a  means  of 
testing  the  techniques  used  for  these  measurements.  In  addition  we  wanted  to  reduce  the  considerable 
uncertainty  associated  with  the  rate  constant  for  this  reaction  (the  reported  values  range  from  0.39x10' 
"  to  3.5x10 11  cm3  s'1). 

A  time-resolved  photolysis  -  probe  technique  is  being  used  to  follow  H  atom  kinetics.  H 
atoms  are  produced  by  pulsed  laser  photolysis  of  H2S.  The  atoms  are  detected  using  two-photon 
excitation  of  the  22S-12S  transition  at  243,2  nm.  Collisions  convert  22S  to  22P,  and  the  subsequent 
Lyman  a  emission  at  121.6  nm  is  detected  by  a  solar  blind  photomultiplier  (Hamamatsu  R6835).  A 
narrow  band  interference  filter  (Acton  122-N-1D)  is  used  to  block  scattered  light  from  the  photolysis 
and  probe  lasers.  The  decay  of  the  H  atom  concentration  is  followed  by  varying  the  delay  between 
the  photolysis  and  probe  pulses. 

Mixtures  of  H2S/F2  or  H2S/C12  flowed  slowly  through  the  photolysis  cell.  The  primary  carrier 
gas  used  in  these  experiments  was  Ar.  To  ensure  that  the  measurements  were  made  under  pseudo  first 
order  conditions  the  H2S  pressure  was  kept  below  0.02  Torr.  F2  or  Cl2  partial  pressures  were  in  the  0.1 
to  0.3  Torr  range.  Ar  buffer  gas  was  added  to  bring  the  total  pressure  in  the  cell  to  around  100  Torr. 
This  was  done  to  effect  translational  cooling  of  the  H  atoms  prior  to  reaction  and  to  limit  diffusion. 
Translational  relaxation  of  the  H  atoms  was  verified  by  examining  the  excitation  profile  of  the  22S-12S 
transition.  The  Doppler  profile  of  this  line  was  consistent  with  a  translational  temperature  of  300  K. 

Kinetic  runs  were  made  with  the  probe  laser  tuned  to  the  center  of  the  excitation  line.  Fig.  10 
shows  a  comparison  of  the  H  atom  decay  following  photolysis  Ar/H2S  and  Ar/H2S/Cl2  mixtures. 
Rapid  removal  of  H  by  the  reaction  with  Cl2  is  easily  observed  by  this  technique.  A  plot  of  the  H  atom 
decay  rate  vs.  [C1J  yielded  a  rate  constant  of  (2.8±0.3)xl0’u  cm3  s'1.  This  value  is  within  the  reported 
range  and  supports  measurements  that  yielded  relatively  large  rate  constants. 
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Figure  10.  Representative  H  +  Cl2  data  is  shown.  The  slow  decay  in  the  absence  of  Cl2  in  the  upper  panel  is 
attributed  to  diffusion  and  H  atom  loss  via  reaction  with  H2S.  The  lower  panel  demonstrates  the  quality  of  the 
fit  to  a  single  exponential  decay.  The  corresponding  rate  constant  is  2.97  ±  0.04  cm3  molecules1  s'1. 
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Attempts  to  characterize  the  H+F2  reaction  by  this  approach  were  less  successful.  There  were 
indications  that  some  of  the  problems  may  stem  from  a  pre-reaction  between  H2S  and  F2.  Previous 
studies  report  that  this  reaction  is  very  slow,  but  adding  F2  to  the  photolysis  cell  dramatically  reduced 
the  signal  from  prompt  H  atoms  (the  signal  seen  at  the  minimum  photolysis-probe  delay  time).  The 
temporal  profile  of  H  atoms  in  the  presence  of  F2  exhibited  multi-exponential  decay  characteristics. 
The  reaction  of  H  with  F2  produces  vibrationally  excited  HF,  and  we  considered  the  possibility  that 
multi-exponential  decay  was  a  consequence  of  the  sequence 

H  +  F2->HF(v)  +  F 
HF(v)  +  hv(243  nm)  — » H  +  F 
H  +  2hv(243)  — > H* 

but  models  based  on  this  scheme  were  unable  to  reproduce  the  complete  family  of  decay  curves.  We 
also  considered  the  possibility  that  a  product  from  the  pre-reaction  was  complicating  the  kinetics 
through  secondary  photochemical  reactions.  To  test  this  notion  we  are  experimenting  with  ways  to 
minimize  the  time  between  mixing  the  reactants  and  photo-initiating  the  reaction.  We  have  added  a 
pulsed  valve  to  the  photolysis  cell  so  that  the  H2S  can  be  introduced  in  brief  busts,  and  the  reaction 
initiated  within  a  few  milliseconds  after  mixing.  Measurements  with  this  system  are  currently  in 
progress. 
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Alexey  L.  Kaledin,  Emory  University,  December  2000 

“Spectroscopy  and  dynamics  of  small  molecules  in  the  gas  phase:  Application  of  theory" 

Interactions 

We  have  kept  in  contact  with  the  group  at  AFRL/DELC  who  are  working  on  the  development 
of  AGIL  (Drs.  Gordon  Hager  and  Gerald  Manke  II).  Our  rate  constants  for  NCl(a)  quenching  and 
self-annihilation  are  being  used  in  their  models  of  the  sub-sonic  AGIL  device,  and  in  the  design  of  a 
super-sonic  device.  The  computational  modeling  group  at  Logicon  RDA  has  identified  key  reactions 
in  the  AGIL  system  for  which  temperature  dependent  rate  constant  and  product  branching  ratio  data 
are  needed.  Where  possible,  we  are  adapting  our  research  program  to  obtain  these  data.  We  remain 
focused  on  the  chemical  kinetic  and  dynamics  issues,  to  ensure  that  our  work  has  intrinsic 
fundamental  value  consistent  with  the  aims  of  the  Molecular  Dynamics  program. 

As  noted  above,  the  kinetic  studies  of  H+C12  and  H+F2  developed  from  conversations  with  the 
group  at  AFRL/DELC.  Dr.  Manke  visited  our  laboratory  for  the  week  of  Sept.  15-21  to  participate  in 
these  experiments.  We  are  currently  preparing  a  manuscript  that  describes  the  H+C12  results. 

The  P.I.  participated  in  a  one-day  workshop  on  the  AGIL  project  at  Logicon  RDA  (NM)  on 
July  23. 

Studies  of  ro-vibrational  energy  transfer  between  the  levels  of  CO(X)  are  being  carried  out  at 
AFRL.  This  project  provides  data  that  may  be  used  to  interpret  and  model  the  behavior  of  CO 
overtone  lasers.  Capt.  Stephen  Phipps  made  an  extensive  study  of  CO(v,J)+CO  collisions  using  IR 
pump  probe  techniques.  As  both  the  excitation  and  probe  lasers  had  sub-Doppler  linewidths,  the  data 
also  revealed  interesting  facets  of  the  competition  between  translational  and  rotational  relaxation  in 
this  system.  Capt.  Phipps  was  registered  with  the  University  of  New  Mexico  (UNM)  as  a  Ph.  D. 
student,  and  the  CO  work  formed  the  core  of  his  graduate  research.  Interpretation  of  the  energy 
transfer  data  was  not  straightforward,  and  Capt.  Phipps  advisors  at  UNM  are  not  experts  in  this  field. 
At  the  request  of  Dr.  Gordon  Hager  I  participated  in  the  interpretation  of  the  data,  and  traveled  to 
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AFRL  on  June  8  for  this  purpose.  I  then  agreed  to  be  a  member  of  the  examining  committee  for  Capt. 
Phipps  Ph.  D.  defense.  I  returned  to  UNM  on  Sept.  17  to  attend  the  thesis  defense,  which  was 
successful.  I  am  continuing  my  involvement  in  this  project,  and  we  are  currently  preparing  a 
manuscript  that  describes  the  key  findings  of  Dr,  Phipps  work. 

We  are  collaborating  with  Dr.  S.  J.  Davis  (Physical  Sciences  Inc.)  on  the  development  of  a 
diode  laser  based  diagnostic  tool  for  characterization  of  NCI(X)  in  AGIL  systems.  Our  most  recent 
contribution  to  this  project  has  been  calculation  of  absolute  linestrengths  for  selected  lines  of  the  b-X 
0-0  band.  These  data  can  be  used  to  convert  quantitative  integrated  absorption  measurements  to 
'  absolute  concentrations.  The  results  of  our  calculations  were  communicated  to  PSI  in  July. 
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